Electron transport in coupled chains of interacting fermions with impurities 
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We study the low-temperature transport of a doped two-chain ladder system of interacting fermions 
in the presence of a barrier or of a low concentration of impurities. Above a certain value of the 
interaction, the conductance is suppressed, like for a single chain, despite the presence of dominant 
superconducting correlations. There is, however, a region of repulsive interaction where perfect 
transmission across the barrier occurs unlike the single-chain case. We provide a possible explanation 
for the temperature maximum of the resistivity in the normal state of Sro.4Cai3.6Cu24 042-5- 



PACS numbers : 71.10.Pm, 72.10.Fk, 71.10.Hf 

Transport properties in one-dimensional (ID) inter- 
acting fermion systems, i.e. Luttinger liquids Q, are 
remarkably affected by the interaction. Indeed, it has 
been shown Q that the linear conductance across an im- 
purity is completely suppressed in the zero-temperature 
limit for repulsive interaction no matter how weak the 
barrier potential. Recently, the question of what hap- 
pens when a small number of Luttinger Liquids are cou- 
pled by a single-particle hopping has received partic- 
ular interest [0-0]. Coupled chains represent the first 
step towards a two-dimensional system, but the impor- 
tant effects of one-dimensionality in the single chains can 
still be treated in an accurate way. Various numerical 
calculations and renormalization-group (RG) treatments 
agree about the fact that the two-chain system is in a 
spin-gapped state with dominant "odd-mode" (the anal- 
ogy of d-wave in two chains) superconducting correla- 
tions |s] The apparent similarity with the high-T c 
superconductors has suggested the idea that the pairing 
mechanism in the high-T c materials could be of a simi- 
lar origin. This theoretical interest is supplemented by 
experimental realizations of systems of 2 to 5 coupled 
chains, such as VO2P2O7, Sr„_iCu„+i02n, and others 
. In addition to the observation of a remnant spin gap 
in the doped metallic state, superconductivity has been 
recently observed at high pressure in the ladder material 
Sro.4Cai3.6Cu24 4 2- ( 5 ||. Finally, narrow GaAs quan- 
tum wires with few transverse channels can also be rep- 
resented by coupled-chain systems in the intermediate 
situation between one and many channels 

In this Letter, we study how the results of Refs. ||] 
for the transport of interacting electrons across an impu- 
rity in a ID Luttinger liquid evolve when two ID sys- 
tems are coupled together. In particular, we analyze 
the influence of the gaps (in the spin and in the "odd" 
charge modes) and of the superconducting correlations, 
present in the two-chain system, on transport proper- 
ties. Specifically, we study the conductance of a system 
of two Hubbard-like chains coupled by a single-particle 
hopping in the presence of a barrier or a low impurity 
concentration. This system is a good description for a 
two-channel conductor and for the ladder materials. Sim- 
ilarly to the pure ID case we obtain that the conduc- 
tance across the impurity is completely suppressed when- 



ever the "even charge" correlation exponent K Pi + is less 
than unity |l0| . However, there is an additional region 
of repulsive electron-electron interaction, U % U c , where 
K p + > 1 and perfect transmission across the impurity 
occurs. A similar result is obtained for the three-chain 
case where, in the pure case, more gapless modes are 
present |]ll| , |r2| . If we limit ourselves to weak impurities 
and temperature T smaller than all the gaps, then our 
results for the T-dependence of the conductivity a agree 
with a recent treatment of a large concentration of weak 
impurities fl3fl . 

We consider a Hamiltonian of the general form 



H — Hchain.t + Hchain.U + H± + Vb 



(1) 



where H cham ^ + H chain .u describes n c (= 2) uncoupled 
ID Hubbard chains, with H c h a in,t denoting the intra- 
chain hopping t, and H c hain,u the on-site repulsion U. 
The chains are coupled by an inter-chain hopping term 
H±_ with hopping amplitude ij_ (with open boundary con- 
ditions in the perpendicular direction) . Vb is the impurity 
potential described in detail below. It should be pointed 
out that the qualitative results do not depend signifi- 
cantly on the detailed structure of the model. One can 
easily extend the method and the results to other systems 
with, e. g., a more general band structure than the simple 
cosine band and a nearest-neighbor interchain and/or in- 
trachain interaction. In the weak-coupling case, one first 
diagonalizes the single-particle part of the Hamiltonian 
(H c hain,t + Hx), obtaining n c bands. The interaction 
Hchain.u written in the band representation includes dif- 
ferent classes of processes that scatter particles or holes 
at the "Fermi surface" (which consists of 2n c points). 
These interactions have to be treated by a RG analysis, 
since perturbation theory is divergent. The clean case 
Vb = 0) has been considered by several authors (Refs. 
[| for two chains and jll|] for three chains). 
One introduces four boson fields ip^ s - q ^> and four dual 
fields 9( s > q *> for the two spins s = ±1 and for the two 
bands q — 0,7r. Taking the symmetric and antisym- 
metric combinations of the band and of the spin part, 
one obtains the fields <p p ,±, <Po± y &p,±, and 9 a _± p0[ . 
The RG calculations of Refs. [|]](| yield a mass gap with 
short-range correlations for both "spin" degrees of free- 
dom <fcr,± and for the "odd-charge" degree of freedom 
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6* p ._ jL4j. In this phase, "odd mode" pairing correlations 
decay rather slowly (as r~ x ^ 2Kp <+^') and dominate with 
respect to other correlations in the weak-coupling limit. 

Next, we introduce the perturbation due to an impu- 
rity in a single chain, or due to a symmetric "barrier" in 
both chains: 



V b = dx 



E 

H,Q2,s 



suppressed due to the presence of the gap in the 9 p - and 
tp a ,- field (at least in phases I and III of Ref. [|1|). Thus, 
they argue that the system should behave as a perfect 
conductor even in the presence of impurities. We show 
here, however, that their result is invalid if one carefully 
takes into account the additional terms in Eq. (Q) gener- 
ated under the RG. Indeed, perturbation terms where the 



Vqi,q2(B)i>t, qi (x)i>s,q 3 ( x ) > ( 2 ) subset of integers n 



where ip Stq (x) destroys a fermion with spin s on band q 
at site x. Here, the impurity potentials V qil q 2 (x) satisfy 
Vo,o(x) = V 7r ^ n (x) and, for a symmetric barrier, Vo,ir{x) = 
0. Introducing the boson representation and considering 
a potential that is restricted to a small region around x = 
0, one obtains three contributions for the perturbation : 
a backscattering intraband term, i.e. 

Vtntra ~ ^ V fi(2k F ) COs((fi p . + + 8<p„ t +) COs((fi p - + Slp a -) 

S 

(3) 

and two interband (back- and forward-scattering) terms 
(the forward-scattering intraband term can, as usual, be 
absorbed by a canonical transformation). The interband 
terms have a form similar to Vintra, except that some 
cosines are replaced with sines and the phase fields 
(for the backscattering) or <p v ,+ (for the forward scatter- 
ing) are replaced with Q v - (v = p or a). In Eq. (||), 
vofl(k) is proportional to the Fourier transform of Vo,o(a;) 
with longitudinal momentum k [we have approximated 
«o,o (2kp) by v g q{2kp{q))], and the fields are implicitly 
considered to be at the origin, x = 0. 

As pointed out in Ref. j| , the perturbation will evolve 
under RG into a more general form consistent with the 
symmetries of the Hamiltonian. In our case, Vb is sym- 
metric under rotation by 2tt of any of the fields ^p v ,+, 
ipv,-, or 9 V< -. Moreover, a simultaneous rotation of ir/2 
of all six of these fields or a simultaneous rotations by 7r 
of some groups of two or three fields (such as, e. g., <p p ,+, 
<p p ,—, and 6cr.-) leaves the perturbation (and also the rest 
of the Hamiltonian) invariant. The general perturbation 
operator which is consistent with these symmetries de- 
velops under the RG and can be cast into the form 



V b '= ]T V(n)exp 



< E 

v=(p,a) 



E n %v»,p 

(4) 



Here, the "Fourier coefficients" V(n) are nonzero only for 
values of the integers n = {njf^} satisfying certain con- 
straints imposed by symmetries. A calculation, in prin- 
ciple similar to ours, has been carried out by T. Kimura 
and coworkers for the two-chain case p5[ . However, these 
authors restrict their analysis to the bare perturbations 
[Eq. (0) plus the interband terms] , which are completely 



and n a _ vanishes survive. The 



"p,+ 



most relevant terms are the ones with the smallest 

consistent with the constraints on the n^f±^ discussed 
above. One (nontrivial) solution is given, for example, 
by the perturbation term v 2 cos(2<^ Pi +) . Here, we have 
exploited inversion symmetry and replaced the gapped 
fields with their average (vanishing) values. To see how 
such a term can originate, consider the intraband part 
of the initial perturbation in Eq. (||) . Although its co- 
efficient vq(E) (whose bare value is vo(Ep) = vo,o{2kp)) 
is suppressed during the RG transformation V2 has 
the RG equation 



dv 2 {E) 
-d \ogE 



l[v (E)] 2 + (l-K p , + )v 2 (E) , (5) 



where E is the characteristic energy parameterizing the 
RG flow. Therefore x>2, although having vanishing ini- 
tial conditions v 2 {E = Ep) — 0, is obtained under 
RG from the irrelevant term vo(E) and is itself rele- 
vant whenever K p _+ < 1. Physically, the fermionic 
representation of the operator associated with v 2 is 

i ; +,s,o( x )" l l ; +,s,7r( x )*P-,s,n( x )*P-,s,o{x) 0, which repre- 
sents an interband electron-pair backscattering operator. 
This term is absent in the bare impurity potential, but is 
obtained to second-order perturbation theory in Vi ntra . 

When K p + > 1, the perturbation v 2 becomes irrel- 
evant and the low-temperature (T smaller than all the 
gaps) conductance G(T) behaves as SG(T) = G(T) 



Go oc -v 2 {Ty 



oc 



_ V 4 T 2(k p ,+-i) m where Go is 



the value of the conductance without impurities, vq = 
vo(Ef), and v 2 (T) is the solution of Eq. (||). For 
K Pt+ < 1, the expression for SG(T) is valid only down 
to a temperature T = T* with SG(T*)/G < 1. For a 
low impurity concentration n, such that the temperature 
is larger than the "discretization energy" T^is = vpn, 
the impurities can effectively be treated independently 
and their contribution to the total resistance adds up 
incoherently 0]. Thus, the resistivity p behaves like 
p-oc v§ n T 2 ( K »-+- 1 ) (down to T = max(T*, T dis )), in 
agreement with Ref. 

When K p .+ < 1 and T < T* , the renormalized pertur- 
bation diverges and one has to treat the case of large v 2 . 
In this limit, and for a single impurity, we can fermion- 
ize back the operator associated with v 2 and consider 
two decoupled semi-infinite (spinless) Luttinger Liquids 
connected via a weak tunneling matrix element t 2 ||. 
This gives a perturbation dual to v 2 : t 2 cos(20 jO , + (O)). 
Similarly to g], one obtains a conductance G(T) oc 
t\ T( 2 / Kp '+)~ 2 , which vanishes at low temperatures. 
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Therefore, the transmission through an impurity in 
a two-chain system is similar to the case of a single 
chain. For K Pj + < 1 the impurity is perfectly insulat- 
ing, whereas perfect transmission occurs in the case of 
K Pt+ > 1, independent of the barrier's strength and sym- 
metry. Physically, this can be understood by considering 
that V2 is associated with a charge-density wave (CDW) 
operator with wave vector 2kp(0) + 2fcp(7r) = Akp whose 
correlations decay like r~ 2K p-+. The argument of Ref. 
M can thus be applied to these 4kp CDW correlations 
(since the "regular" 2Uf ones are suppressed jl5| ) . When 
K p , + < 1 the Akp CDW correlations are long-range 
enough and are pinned by an arbitrarily weak impurity 
at low energies. 

Nevertheless, there is an essential difference between 
the single- and the two-chain system, since if Pl + is renor- 
malized by an amount SK Pt+ oc (Svf/vf) 2 by the pres- 
ence of the gap A Q in the "odd" modes || (here, Svf 
is the difference between the Fermi velocities dj? of the 
bonding and antibonding bands). This renormalization 
can be considered as an effective attraction mediated by 
the odd gapped modes. The RG equation for K Pt+ as 
a function of the flux parameter T (T is proportional 
to I = — XogE/Ep in our notation) is reported in Ref. 
Q (in order to match our notation, one has to replace 
K p + — > l/K Pt+ in that reference). The exact value of 
8K Pi+ cannot be determined even in the U/t — > limit, 
since it depends on the point I* where the RG flow is 
stopped. A reasonable lower bound and estimate for 
SK Pj+ and for the critical value U c , which separates re- 
gions K p _ + < 1 and K p>+ > 1 in Fig. 1 can be obtained in 
the following way: Notice that 8K p + contains two con- 
tributions: (i) an /-independent contribution obtained by 
neglecting the RG flow of the couplings (i. e. by taking 
/ = and by integrating out the gapped modes), and (ii) 
an additional contribution due to the renormalization of 
the couplings which depends on I. One can show that (i) 
sets a lower bound for SK Pi+ and that the contribution of 
(ii) is much smaller than (i) and only weakly /-dependent. 

In Fig. 1, we plot U = U c as a function of doping S 
(away from half-filling) for the isotropic case t± = 1, by 
taking into account in the dashed line only the contribu- 
tion (i) and by including in the full line additionally the 
contribution (ii) jl7) . For large doping there is a sizable 
region of repulsive interaction for which the transmission 
across the impurity is perfect. This effect will be particu- 
larly relevant for systems in which the Hubbard repulsion 
U < t, as for example in carbon nanotubes Jl8| and in 
GaAs quantum wires. Exact numerical results for small 
U 1 19 yield the opposite effect, namely K p _ + seems to be 
reduced with respect to the single-chain value. We be- 
lieve, however, that in these numerical studies the value 
of K Pj+ is underestimated since the system size consid- 
ered (6 lattice spacings oq) is smaller than the spin cor- 
relation length £ s oc Uf/A s , which gives £ s > 30 — 40ao 
P], where A s is the spin gap. 

For a finite impurity concentration, the resistivity p(T) 
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FIG. 1. Phase diagram for the two-chain system with a 
barrier or with a low impurity concentration. The two esti- 
mates for the boundary U — U c are obtained according to 
contribution (i) only (dashed) and additionally (ii) (full) as 
explained in the text. 



diverges at low temperatures in the region in which 
K p + < 1, whereas it decreases in the region in which 
K p , + > 1 pb|| . At temperatures larger than the gap, 
T 3> A Q w A s p]]| , the resistivity is dominated by the 
scattering with the 2fcp potential vq. Its behavior (2^] is 
p(T) oc v 2 jiT <yK ''~ 1 ^l 2 , yielding an increasing resistivity 
with decreasing T (for U > 0) . 

For a more quantitative evaluation of p(T) and in or- 
der to find the behavior of p(T) in the crossover regime 
for T of the order of A„, we use the perturbative ex- 
pression p(T) oc c v (T) 2 + C2 V2(T) 2 , where cq and 
C2 are constants J23|. For these temperatures, v (T) 
can be described by an activated behavior vq(T) « 
(T/^ F )(^- 1 )/ 4 u exp(A o / J B F - A /T), whereas v 2 can 
be evaluated from Eq. (||) p2| . In the weak-impurity 
regime which we are considering, the parameter v = 
vo \/c2 / Co Q takes a value v < 1. In this case p(T) 
is dominated by the 2fc F -scattering potential vq from 
higher temperatures down to a crossover temperature 
T = T m < A , with T m - A /(21og(2/«)). Therefore, 
p(T) eventually decreases down to this crossover temper- 
ature T m . For T < T m the 4fc F -scattering v 2 becomes 
dominant and, if K p , + < 1, p(T) dramatically increases 
eventually diverging for low T. 

In Fig. 2, we plot the behavior of p{T) for differ- 
ent values of the parameters v and A S /A Q , and for 
K p w K p . + = 0.6, which corresponds to U/t » 8.0 |J] 
and doping S ~ 0.25 [|o|j25|. The full line is obtained 
by taking a large ratio A S /A D w 4 between the spin 
and the odd-mode gaps. For this large ratio and in the 
regime 2A Q < T < 2A S , p(T) is dominated by the factor 
oc T( Kp ~ 3 ^ 2 , since in this regime the fields fa,± have a 
vanishing expectation value |2(| . This induces the down 
turn of p(T) at these higher temperatures. 

The results shown by the full line in Fig. 2 may 
provide one possible explanation of the experimentally 
observed maximum of p(T) at T « 100 K in the nor- 
mal state of Sro.4Cai3. 6 Cu24042_(5 [§,^4|- As mentioned 
above the large ratio A S /A Q guarantees the experimen- 
tally observed down turn of p(T) above its maximum. 
Furthermore, our choice of A S /A Q w 4 is consistent 
with the following estimate: from the full line in Fig. 
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2, A s has a value of about 1.3 times the temperature 
for which p(T) has a maximum. Comparing with ex- 
periment, this fixes A s , i.e. A s « 130 K. On the other 
hand, if we take the value for the undoped material, 
A s (n — 1) « 430 K ||], and estimate the doping-induced 
reduction (according to the second of Ref. we ob- 
tain A s (n = 0.75) ks 0.33A s (?i = 1) » 140 K. This 
is consistent with our previous result. The increase of 
p(T) at T > 200 K in the experiments is probably due to 
contributions from electron-electron and electron-phonon 
scatterings. Notice that an appropriate description of 
the superconducting transition cannot be obtained by our 
simplified effectively one-dimensional model. We finally 
comment on the pressure dependence of p(T) shown by 
the experiments. Application of hydrostatic pressure in- 
creases the carrier density in the ladders p7j| , thus de- 
creasing the overall resistivity (as observed experimen- 
tally), and at the same time slightly reducing the values 
of the gaps and, consequently, the peak position in p(T) . 
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4.0 . . 6.0 

T/A 

FIG. 2. Behavior of the resistivity p as a function of 
temperature (in units of the odd-mode gap A„) for values 
of the parameters (as explained in the text): v = 0.4 and 
A s /A = 4. (full line); v = 0.4 and A s = A D (dotted); v = 0.1 
and A s = A D (dashed). For all curves, K p = K Pt+ = 0.6 
and _Ef/A d = 150. Inset: p{T) for the normal state of 
Sri4_ :E Ca :E Cu24 42-5 from Ref. [8] for different values of pres- 
sure. (With kind permission of the authors). 

In conclusion, we have studied the transport properties 
of two-chain ladder systems in the presence of a barrier or 
of a low impurity concentration. We find a small region 
of repulsive interaction < U < U c [cf . Fig. [jj , for which 
the system is a perfect conductor for Tdi S ^ T < A G . We 
also suggested a possible explanation for the normal state 
behavior of the resistivity in Sri4_ x Ca a; Cu24042-5 . 
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